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Abstract. We calculate static properties of non- rotating 
neutron stars (NS's) using a microscopic equation of state 
(EOS) for asymmetric nuclear matter, derived from the 
Brueckner-Bethe-Goldstone many-body theory with ex- 
plicit three-body forces. We use the Argonne AV14 and 
the Paris two-body nuclear force, implemented by the Ur- 
bana model for the three-body force. We obtain a maxi- 
mum mass configuration with Mmax = 1.8Mq {Mmax = 
I.94M0) when the AV14 (Paris) interaction is used. They 
are both consistent with the observed range of NS masses. 
The onset of direct Urea processes occurs at densities 
n > 0.65 fm^^ for the AV14 potential and n > 0.54 fm~^ 
for the Paris potential. Therefore, NS's with masses above 
j^Urca ^ 1 for the AV14 and M^'''"' = 1.2AMq for 

the Paris potential can undergo very rapid cooling, de- 
pending on the strength of superfluidity in the interior of 
the NS. The comparison with other microscopic models 
for the EOS shows noticeable differences. 

Key words: three-body forces-equation of state-neutron 
stars 



1. Introduction 

In the next few years it is expected that a large amount of 
novel informations on neutron stars (NS's) will be avail- 
able from the new generation of X-ray and 7-ray satel- 
lites. Therefore, a great interest is devoted presently to 
the study of NS's and to the prediction of their structure 
on the basis of the properties of dense matter. The equa- 
tion of state (EOS) of NS matter covers a wide density 
range, from ~ 10 g/cm^ in the surface to several times 
nuclear matter saturation density (po 2.8 10"^^ g/cur') 
in the center of the star (Shapiro et al. 1983). The interior 
part (core) of a NS is made by asymmetric nuclear matter 
with a certain lepton fraction. At ultra-high density, mat- 
ter might suffer a transition to other exotic hadronic com- 
ponents (like hyperons, a K~ condensate or a deconfined 



phase of quark matter). The possible appearance of such 
an exotic core has enormous consequences for the neutron 
star and black hole formation mechanism (Bombaci 1996). 
Unfortunately large uncertainities are still present in the 
theoretical treatment of this ultra-dense regime (Glenden- 
ning 1985, Prakash et al. 1997). Therefore, in the present 
work, we consider a more conventional picture assuming 
the NS core is composed only by an uncharged mixture 
of neutrons, protons, electrons and muons in equilibrium 
with respect to the weak interaction (/3-stable matter). 
Even in this picture, the determination of the EOS of 
asymmetric nuclear matter to describe the core of the NS, 
remains a formidable theoretical problem (Hjorth- Jensen 
et al. 1995). 

Any "realistic" EOS must satisfy several requirements 
: i) It must display the correct saturation point for sym- 
metric nuclear matter (SNM); ii) it must give a symmetry 
energy compatible with nuclear phenomenology and well 
behaved at high densities; iii) for SNM the incompressibil- 
ity at saturation must be compatible with the values ex- 
tracted from phenomenology (Myers et al. 1996); iv) both 
for neutron matter (NEM) and SNM the speed of sound 
must not exceed the speed of light (causality condition), 
at least up to the relevant densities; the latter condition 
is automatically satisfied only in fully relativistic theory. 

In this work we present results for some NS properties 
obtained on the basis of a microscopic EOS which satisfies 
requirements i-iv, and compare them with the predictions 
of other microscopic EOS's. 



2. Equation of state 

2.1. Brueckner-Bethe-Goldstone theory 

The Brueckner-Bethe-Goldstone (BBG) theory is based 
on a linked cluster expansion of the energy per nucleon of 
nuclear matter (see e.g. Bethe 1971). The basic ingredi- 
ent in this many-body approach is the Brueckner reaction 
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matrix G, which is the solution of the Bethe-Goldstone 
equation 

\kakb)Q{kakb\ 



G{n;Lo)=v + vJ2^ 



kakb 



w - e{ka) - e{kb) 



G{n;u}), 



(1) 



where v is the bare nucleon-nucleon (NN) interaction, n 
is the nucleon number density, and u the starting energy. 
The operator \kakb)Q{kakb \ projects on intermediate scat- 
tering states in which both nucleons are above the Fermi 
sea (Pauli operator). 



e(fc) = e(fc; n) 



*2 

— k'^ + U(k;n) 
2m 



(2) 



is the single particle energy. The Brueckner-Hartree-Fock 
(BHF) approximation for the s.p. potential U {k; n) using 

the continuous choice is 

U{k;n)= ^ {kk'\G{n;e{k) + e{k'))\kk')a (3) 

k'<kF 

where the subscript "a" indicates antisymmctrization of 
the matrix element. In this approach equations (l)-(3) 
have to be solved selfconsistently. In the BHF approxima- 
tion the energy per nucleon is: 
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(A:fc'|G(n; e(fc) + e{k'))\kk')a (5) 



k,k'<kF 



In this scheme, the only input quantity we need is the 

bare NN interaction v in the Bcthc-Goldstonc equation 
(1). In this sense the BBG approach can be considered as 
a microscopic one. 

The Brucckncr-Hartrcc-Fock (BHF) approximation for 
the EOS in SNM, within the continuous choice (Baldo et 
al. 1991), reproduces closely results which include up to 
four hole line diagram contributions to the BBG expansion 
of the energy, calculated within the so called gap choice for 
the single particle potential U{k) (Day et al. 1985). In the 
following we refer to the latter calculations as BBG^'^^. 
The numerical accuracy of the solution of the BHF equa- 
tion in the continuous choice has been discussed recently 
in ref. (Schulzc et al. 1995). 

In Fig. 1 we show the energy per nucleon calculated 
within this scheme in the case of SNM (lower curves) and 
NEM (upper curves); the solid and dashed curves have 
been obtained using respectively the Paris potential (La- 
combe et al. 1980) and the Argonne Vi4 (AV14) model 
(Wiringa ct al. 1984) for the two-body nuclear force. For 
comparison, we report in the same figure the energy per 
nucleon calculated within the variational many-body ap- 
proach (Wiringa ct al. 1988, hereafter WFF) using re- 
spectively the AV14 (open circles) and the Urbana UV14 
(squares) two-nucleon potentials. 

We notice that both the BHF and the variational ap- 
proach fail to reproduce the empirical saturation point of 
nuclear matter (n^ = 0.17 ±0.01 fm-^, Eo/A = -16 ± 
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Fig. 1. The energy per baryon E/A is plotted vs. the num- 
ber density n for symmetric matter (lower curves) and for 
neutron matter (upper curves). The solid (dashed) line 
represents a non-relativistic Brueckner calculation with 
Paris (AV14) potential, whereas the open circles (squares) 
are the results of a variational calculation performed with 
AV14 (UV14) interaction. 



1 MeV). Analogous non-relativistic calculations have been 
performed with the Bonn potential (Machleidt 1989, En- 
gvik et al. 1996). Even in that case symmetric nuclear 
matter does not saturate at the empirical point. 

The comparison of the saturation curves obtained in 
the BHF and the variational approach with the same two- 
nucleon AV14 potential shows that our BHF saturation 
point (no = 0.256 /m'^, Eo/A = -18.26 MeV) is closer 
to the empirical saturation point than the one reported 
in WFF paper {uo = 0.319 /W^, Eo/A = -15.6 MeV). 
The two methods give results in reasonable agreement in 
the case of SNM up to a density n ~ 0.4 ,fm~^, whereas 
sizeable differences are evident for NEM already at densi- 
ties n ~ 0.3 fm~^. A similar outcome has been found by 
comparing variational and BBGf^^ saturation curves for 
the same NN potential (Day, et al 1985). The discrepancy 
clearly arises from the different many-body technique. 

In Fig. 2 we report the energy per nucleon as a func- 
tion of the number density n calculated in the BHF ap- 
proximation using the Paris (solid curves) and the AV14 
potentials (dashed curves) up to densities typically en- 
countered in the core of a neutron star. The upper curves 
refer to neutron matter and the lower ones to symmetric 
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nuclear matter. We notice that the EOS derived with the 
Paris potential is stiffer than the one obtained with the 
AV14 interaction. This is possibly due to the momentum 
dependence of the Paris potential, as discussed recently in 
ref.(Engvik et al. 1997). 
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Fig. 2. The energy per baryon E/A is plotted vs. the num- 
ber density n for symmetric matter (lower curves) and for 
neutron matter (upper curves). Non-relativistic Brueckner 
calculations are drawn for the AV14 (dashed line) and the 
Paris (solid line) potentials. 



2.2. Three-body forces 



Non-relativistic calculations, based on purely two-body 
interactions, fail to reproduce the correct saturation point 
for symmetric nuclear matter (Coester et al. 1970). This 
well known deficiency is commonly corrected introducing 
three-body forces (TBF). Unfortunately, it seems not pos- 
sible to reproduce the experimental binding energies of 
light nuclei and the correct saturation point accurately 
with one simple set of TBF (Wiringa et al. 1988). Relevant 
progress has been made in the theory of nucleon TBF, but 
a complete theory is not yet available. A realistic model 
for nuclear TBF has been introduced by the Urbana group 
(Carlson et al. 1983; Schiavilla et al. 1986). The Urbana 
model consists of an attractive term V^^ due to two-pion 



exchange with excitation of an intermediate A-resonance, 
and a repulsive phenomenological central term V^^^ 

(6) 



ijk 



^ijk 



The two-pion exchange contribution is a cyclic sum over 
the nucleon indices i, j, k of products of anticommutator 
{,} and commutator [,] terms 



eye 



(7) 
(8) 

(9) 

is the one-pion exchange operator, a and r are the Pauli 

spin and isospin operators, and Sij = 3[((Ti ■rij){aj -rij) — 
(Tiaj\ is the tensor operator. Y{r) and T(r) are the Yukawa 
and tensor functions, respectively, associated to the one- 
pion exchange, as in the two body potential. 
The repulsive part is taken as 



where 

Xij = Y{rij)ai ■ aj + T{rij)Sij 



ijk 



U 



eye 



T\ni)T\rjk) 



(10) 



The constants A and U in the previous equations can be 
adjusted to reproduce observed nuclear properties. Within 
this scheme Schiavilla et al. (1986) found A = —0.0333 and 
U = 0.0038 by fitting properties of light nuclei ('^H, '*He), 
the so called Urbana VII parametrization (UVll) (Carlson 
et al. 1983). 

We introduced the same Urbana three-nucleon model 
within the BHF approach. To incorporate the TBF in the 
Brueckner scheme we followed the method of Lejeune et 
al. (1986). The TBF is reduced to an effective two-body 
force by averaging on the position of the third particle, 
assuming that the probability of having two particles at a 
given distance is reduced according to the two-body cor- 
relation function. The resulting effective two-body force is 
of course density dependent. Further details will be given 
elsewhere (Baldo et al. in preparation). We have adjusted 
the parameters A and U in order to reproduce closely 
the correct saturation point of SNM, since for NS stud- 
ies this is an essential requirement, and there is no reason 
to believe that TBF be the same as in light nuclei. Our 
values for the TBF parameters are A' = —0.0329 and 
U' = 0.00361 both for the Paris and the AV14 potentials. 

The corresponding EOS's obtained using the Paris (so- 
lid line) and the AV14 potentials (short dashed line) are 
depicted in Fig. 3 for symmetric (panel (a)) and neutron 
matter (panel (b)). The corresponding values are reported 
in Table 1. Those EOS's saturate respectively at 

Ho = 0.176 fm,-^, Eo/A = -16.01 MeV Paris (11) 

and 

no = 0.178 fm-^, Eo/A = -16.46 MeV AVU (12) 
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Fig. 3. The energy per baryon E/A is plotted vs. the number density n for symmetric matter (panel (a)) and for 
neutron matter (panel (b)). Several EOS's are shown, i.e. non-relativistic Brueckner calculations with three-body 
forces (solid line, Paris potential; short dashed line, AV14 potential) and a relativistic Dirac-Brueckner one (long 
dashes, DBHF). For comparison a variational calculation (with AV14-I-UVII interaction) is also reported (open circles, 
WFF). 



They are characterized by an incompressibility K^"^ - 
281 MeV and K^'^^ = 253 MeV, the latter being very 
close to the recent phenomenological estimate of Myers 
(1996). In the same figure, we show the EOS obtained 
within the variational many-body approach by WFF when 
the UVII parametrization (Carlson et al. 1983, Schiavilla 
et al. 1986) is added to the AV14 two-body force (open 
circles). 

We have already noticed that both the BHF and the 
variational many-body methods produce a different sat- 
uration point already at the two-body level. As a con- 
sequence, the values of the strengths {A' and U') of the 
attractive and repulsive TBF needed to reproduce empir- 
ical saturation are different with respect to those used by 
WFF in their work. In particular, being U' < U, our re- 
pulsive TBF is weaker. Consequently, our EOS's are softer 
than the WFF at very high density (see Fig. 3) where the 
repulsive component of the TBF is dominant. 

In Fig. 3 we plot also the EOS from a recent Dirac- 
Brueckner calculation (DBHF) (Li et al. 1992) with the 
Bonn A two-body force (long dashed line) . In the low den- 
sity region (n < 0.4 fm~^), both BHF equations of state 
with TBF and DBHF equations of state are very similar, 
whereas at higher density the DBHF is stiffcr. The discrep- 
ancy between the non-relativistic and relativistic calcula- 
tion of the EOS can be easily understood by noticing that 
the DBHF treatment is equivalent (Baldo et al. 1995) to 
introduce in the non-relativistic BHF the three-body force 
corresponding to the excitation of a nucleon-antinucleon 



pair, the so-called Z-diagram (Brown et al. 1987). The lat- 
ter is repulsive at all densities. In BHF treatment, on the 
contrary, both attractive and repulsive three-body forces 
are introduced, and therefore a softer EOS can be ex- 
pected. It should be noticed that, because of the strong 
repulsion at short distances, our BHF calculation with 
three-body force is numerically stable up to densities n ~ 
0.8 fm~^. Moreover, the relativistic DBHF calculations 
are available in literature up to densities n = 0.75 fm~^ 
(Li et al., 1992). Since the typical densities encountered in 
the core of neutron stars are larger, both microscopic cal- 
culations need to be either extrapolated to slightly higher 
densities or joined with other high density equations of 
state, several of which are available in literature. We de- 
cided to perform a numerical extrapolation to the needed 
higher densities. In the extrapolation procedure we keep 
checking that the causality condition is fulfilled, i.e. 

c./c.(-) <1, (13) 

which means that the speed of sound Cg in matter must 
be lower than the speed of light in the vacuum. This is a 
basic trait of any "realistic" EOS, regardless the details of 
the interactions among matter constituents or the many- 
body approach. Details on the extrapolation procedure are 
given in Sec. 3. 

The ratio c^/c as a function of the number density is 
reported in Fig. 4, in the case of pure neutron matter, for 
our BHF models (AV14 -|- TBF(short dashes), Paris + 
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Table 1. The energy per baryon E/A in MeV is shown vs. the number density n in fm~^ for symmetric matter and 
neutron matter. The values are the results of a non-relativistic BHF calculation with AV14 and Paris potentials with 
three-body forces. 



n 


AVU + TBFsym 


Paris + TBF^ym 


AVlA + TBFr,^^ 


Paris + TBF„eu 


0.08 


-11.47 


-10.8 


8.54 


9.47 


0.16 


-16.32 


-15.88 


16.2 


17. 


0.2 


-16.25 


-15.74 


21. 


22.1 


0.3 


-10.78 


-9.4 


35.6 


38.8 


0.4 


0.79 


3.2 


55.26 


61.4 


0.5 


17.53 


20.46 


81.13 


90.3 


0.6 


38.5 


41.78 


112.5 


125.55 


0.7 


62.7 


67.17 


147.9 


166.7 


0.8 


89.3 


96.7 


185.9 


212.66 


0.9 


117.67 


130.22 


225.7 


262.3 


1. 


147.35 


167.4 


266.7 


314.6 


1.1 


178. 


207.76 


308.4 


368.78 


1.2 


209.5 


250.87 


350.8 


424.12 


1.3 


241.58 


296.27 


393.5 


480.3 


1.4 


274.17 


343.6 


436.6 


537. 


1.5 


307.18 


392.52 


479.8 


594.2 



TBF (soUd)), and the DBHF model (long dashes). WFF 
model (open circles) violates the causality condition at 
densities encountered in the core of neutron stars near 
the maximum mass configuration predicted by that model 
(Wiringa et al.,1988). 




0.5 



n (fm--^) 



Fig. 4. The ratio Cg/c is plotted as a function of the num- 
ber density for pure neutron matter. The notation is the 
same as in Fig.3. 



2.3. Symmetry energy and EOS for P -stable matter 

The energy per nucleon of asymmetric nuclear matter 
(n„ ^ Hp, being r?,„ and Up the neutron and proton num- 
ber densities respectively) can be calculated extending the 
original BEG theory for SNM sketched in sect. 2.1. (Bom- 
baci et al. 1994). From the energy per nucleon one defines 
the nuclear symmetry energy as 



_ 1 d^E/A 



2 dp^ 



/3=0 



where we introduce the asymmetry parameter 



Tlfi Tip 

n 



(14) 



(15) 



with n = Un+rip. At the saturation density, Egyj^ij^o) cor- 
responds to the symmetry coefficient in the Weizsacker- 
Bethe semiempirical mass formula, and its value extracted 
from nuclear systematics lies in the range 28-32 MeV. 

In the present work, we calculate the energy per nu- 
cleon of asymmetric nuclear matter assuming the so called 
parabolic approximation 



^ (ri, 0) = ^{n,p = 0)+ Esym{n)0' 



(16) 



A' ' A' 

In this approximation the symmetry energy can be ex- 
pressed in tc;rms of the difference of the energy per particle 
between neutron (/3 = 1) and symmetric (/? = 0) matter: 



E. 



sym 



(n) = |(n,/3=l)-|(n,/? = 0) 



(17) 



The reliability of the parabolic approximation has been 

checked microscopically by Bombaci and Lombardo (1991). 
Then we determine the composition of /3-stable matter, 
at a given nucleon number density, solving the equations 
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Fig. 5. The symmetry energy and the proton fraction are shown vs. number density respectively in panel (a) and (b). 
The notation is the same as in Fig. 3. 



for chemical equilibrum (for neutrino-free matter, = 
M.. = 0) 

Me = - i"p = A (18) 

Me = (19) 

and charge neutrality 

«p = "e + n^, (20) 

being //^ {i = n,p, e, /i) the chemical potentials of star con- 
stituents. Therefore, in the case of neutrino-free matter, 
the leptonic chemical potential is determined once the dif- 
ference /i = /i„ — /Xp is known. The latter quantity can be 
expressed as: 



M = 



dE/A 



dx 



= 2 



dE/A 



dp 



(21) 



where x = np/n — (1 — j3)/2 is the proton fraction and the 
partial derivatives are taken for constant nucleon number 
density n. In eq. (21) we neglect the neutron-proton mass 
difference. 

In the parabolic approximation (16) for the energy per 
nucleon of asymmetric nuclear matter one has 

^^AEsy,n{n){l^2x) (22) 

Therefore the composition of /3-stable matter, and in par- 
ticular the proton fraction a;, is strongly dependent on the 
nuclear symmetry energy. 

As it has been recently pointed out by Lattimer et al. 
(1991), the value of the proton fraction in the core of NS 



is crucial for the onset of direct Urea processes, whose 
occurrence enhances neutron star cooling rates. 

The values of Esym for the different EOS's are re- 
ported in Fig. 5 (panel (a)), together with the correspond- 
ing proton fraction (panel (b)). The long dashed line rep- 
resents the relativistic DBHF calculation^, whereas the 
non-relativistic Brueckner calculations with the AV14 4- 
TBF and the Paris -I- TBF models are respectively repre- 
sented by a short dashed and a solid line. We notice that 
in both relativistic and non-relativistic Brueckner-type 
calculations, the proton fraction x can exceed the "criti- 
cal" value x'^''™ = (11 — 15)% needed for the occurrence 
of direct Urea processes (Lattimer et al. 1991). This is at 
variance with the WFF variational calculation (open cir- 
cles), which predicts a low absolute value both for the sym- 
metry energy and the proton fraction with a slight bend 
over. For the AV14 + TBF model we find x'^'^™ = 13.9%, 
which correspond to a critical density n^^'^°- = 0.65 fm~^, 
whereas the Paris + TBF model predicts the same criti- 
cal proton fraction at density n'^''™ = 0.54 fm~^ . There- 
fore, BHF neutron stars with a central density higher than 



^ In a recent paper (Lee et al., 1997) the authors show micro- 
scopic calculations of asymmetric nuclear matter in the DBHF 
approach up to densities n = 1.1 fm~^. They obtain a symme- 
try energy which saturates with increasing density, at variance 
with our findings. In spite of that the values of the maximum 
mass configuration are similar to those calculated in this work. 



M. Baldo et al.: Microscopic nuclear equation of state with three-body forces and neutron star structure 



7 



600 



> 



< 



400 



200 - 




n (fm ) 



Fig. 6. The energy per baryon E/A of asymmetric nu- 
clear matter in /3-equilibrium (with electrons and muons) 
is plotted vs. the number density n. Different curves are 
denoted as in Fig. 3. 

^Urca develop inner cores in which direct Urea processes 
are allowed. 

In Fig. 6, we show the energy per baryon of asymmet- 
ric nuclear matter in /3-equilibrium, for the models under 
consideration. Wo notice^ that the AV14 and the Paris two- 
nucleon interactions, implemented by the same three-body 
force, produce similar EOS's up to densities n ~ 0.8 fm^^ , 
whereas the high density region is quite different. There- 
fore we find the same trend as observed in Fig. 2. 

From the energy per baryon of asymmetric nuclear 
matter in /3-equilibrium, we calculate the nuclear contri- 
bution Pnuci to the total pressure of stellar matter using 
the thermodynamical relation 

^2d{E/A) 



Pnucl (^) 



dn 



(23) 



Then the total pressure and total mass density p are given 

by 

P = Pnucl + Plep (24) 

P = ^ (^nucl + £(ep) = ^ + " + ^'ep) (25) 

being Piep and eiep the leptonic contributions to the total 
pressure and energy density, rriN the nucleon mass and c 
the speed of light in the vacuum. 

Our EOS's for /3-stable matter (neutrons, protons, elec- 
trons and muons) calculated with AV14 (Paris) potential 



Table 2. EOS for /J-stable matter obtained in BHF ap- 
proximation using the Argonne AV14 two-body interac- 
tion implemented by TBF as described in the text. We 
display the number density n in units of fm~^, the proton 
fraction x, the total mass density p and the total pressure 
P. 



n 


X 


p (10'* g cm-") 


P (10^* dyn cm-'') 


0.08 


0.025 


1.35 


0.042 


0.16 


0.048 


2.7 


0.244 


0.2 


0.06 


3.4 


0.502 


0.3 


0.078 


5.15 


1.88 


0.4 


0.093 


6.98 


4.62 


0.5 


0.11 


8.92 


8.92 


0.6 


0.129 


10.97 


14.7 


0.7 


0.148 


13.13 


21.87 


0.8 


0.166 


15.42 


30.36 


0.9 


0.182 


17.85 


40.14 


1. 


0.195 


20.39 


51.19 


1.1 


0.207 


23.07 


63.5 


1.2 


0.218 


25.88 


77.04 


1.3 


0.227 


28.82 


91.8 


1.4 


0.236 


31.88 


107.78 


1.5 


0.243 


35.07 


124.96 


Table 3. Same 


as Table 2, but for Paris potential a 


interaction. 






n 


X 


p (10'* g cm-') 


P iW" dyn cm~'^) 


0.08 


0.025 


1.35 


0.039 


0.16 


0.05 


2.7 


0.243 


0.2 


0.062 


3.4 


0.53 


0.3 


0.085 


5.17 


2.07 


0.4 


0.105 


7. 


4.83 


0.5 


0.128 


8.95 


8.95 


0.6 


0.155 


11.0 


14.8 


0.7 


0.18 


13.2 


23. 


0.8 


0.204 


15.55 


34. 


0.9 


0.223 


18.06 


48. 


1. 


0.237 


20.77 


64.9 


1.1 


0.25 


23.67 


84.53 


1.2 


0.257 


26.78 


106.6 


1.3 


0.264 


30.1 


131. 


1.4 


0.269 


33.6 


157.6 


1.5 


0.272 


37.3 


186.4 



plus Urbana three-body forces are given respectively in 
Table 2 (3). 



3. Neutron star structure 

The EOS for /3-stable matter can be used in the Tolman- 
Oppenheimer Volkoff (Tolman 1934, Oppenheimer et al. 
1939) equations to compute the neutron star mass and 
radius as a function of the central density. 

We assume that a neutron star is a spherically sym- 
metric distribution of mass in hydrostatic equilibrium. 
We neglect the effects of rotations and magnetic fields. 
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Then the equihbrium configurations are simply obtained 
by solving the Tolman-Oppenheimer-Volkoff (TOV) equa- 
tions for the total pressure P and the enclosed mass m, 



dP{r) _ Gm{r)p{r) (1 + s^(fy)(l + ""cMr^ ) 



dr 



dm{r) 

dr 



(1 



2Gm(r) ^ 
rc^ ^ 



= 4:TTr'^p{r) 



(26) 



(27) 



being G the gravitational constant. Starting with a central 
mass density p{r = 0) = pc, we integrate out until the 
pressure on the surface equals the one corresponding to 
the density of iron. This gives the stellar radius R and the 
gravitational mass is then 

Mg = m(i?) =47r / drr^p{r). (28) 



For the outer part of the neutron star we have used the 
equations of state by Feynman-Metropolis- Teller (Feyn- 
man ct al. 1949) and Baym-Pcthick-Suthcrland (Bayni ct 
al. 1971), and for the middle-density regime (0.001 fm~^ < 
n < 0.08 fm~^) we use the results of Negele and Vau- 
therin (Negclc ct al. 1973). In the high-density part (n > 
0.08 fm~^) we use alternatively the four EOS's discussed 
above. The results are reported in Fig. 7. We display the 
gravitational mass Mg, in units of the solar mass Mq 
(Mq = 1.99 10^'^ g), as a function of the radius R (panel 
(a)) and the central number density ric (panel (b)). The 
notation is the same as in the previous figures. 

As expected, the stiffest EOS within the Brueckner 
scheme (DBHF) we used in the present calculation gives 
higher maximum mass and lower central density with re- 
spect to the non-relativistic Brueckner models. 

The properties of the maximum mass configuration, 
for each EOS, are summarized in table 4, whereas in ta- 
ble 5 the properties of neutron stars with gravitational 
mass Mq = 1.4 Mq are reported. However, the proper- 
tics of each configuration extracted by solving the TOV 
equations are slightly dependent on the extrapolation pa- 
rameters we used. This holds for the DBHF equation of 
state and the Paris + TBF model, while in the case of 
Argonne AV14 -|- TBF model no seizable dependence was 
found. For the BHF calculation with only two-body forces 
no extrapolation is needed. In any case the dependence on 
the extrapolation procedure turns out to be rather weak. 
We have used a fractional polynomial extrapolation of the 
form 



E 



g-(n) 



l + 6n^ 

where the coefficients of the polynomial Q are fitted to 

reproduce the EOS up to the maximum calculated den- 
sity n, and the coefficient b ensures that for very large 
values of n the energy per baryon does not diverge. Dif- 
ferent values for the polynomial degree a can be used and 
different extrapolated EOS are in general obtained. How- 
ever we found that the critical mass, radius and central 



Table 4. Properties of the maximum mass configuration 
obtained for different equations of state: Mq is the gravi- 
tational (maximum) mass, R is the corresponding radius, 
ric the central number density and Xc the central proton 
fraction. 



EOS 


Mg/Mq 


R(km) 




) 




AV14 


1.5 


8.1 


1.9 




0.17 


Paris 


1.66 


8. 


1.88 




0.26 


AV14 + TBF 


1.8 


9.7 


1.34 




0.23 


Paris + TBF 


1.94 


9.54 


1.33 




0.265 


DBHF 


2.1 


10.6 


1.07 




0.24 


WFF 


2.130 


9.40 


1.25 




0.045 


Table 5. Properties of neutron stars with Mg 


= 1.4 A 


EOS 


R(km) 


nc(/m-^) 


Xc 






AV14 


8.92 


1.27 


0.12 






Paris 


9.3 


1.08 


0.178 






AV14 + TBF 


11. 


0.65 


0.139 






Paris -1- TBF 


11.09 


0.64 


0.165 






DBHF 


11.74 


0.49 


0.14 






WFF 


10.41 


0.66 


0.0658 







■0 



densities do not vary more than 5% with cr = 3, 4 and b 
ranging from lO^'^ to 10^"^ for the case of the Paris po- 
tential -I- TBF model calculations, while in all other cases 
the variations can be considered ncgligcablc. 

As one can sec from Tabic 4, the values of the limiting 
mass configiiration arc strongly modified by the presence 
of the nuclear TBF in the equation of state, in particular 
by their repulsive contribution which dominates at high 
densities. However, the non-relativistic BHF calculations 
implemented by three-body force do not strongly depend 
on the two-body interaction, and predict very close lim- 
iting values of the maximum mass (within the numerical 
uncertainties due to the extrapolation). 

The difference between BHF and WFF neutron stars 
reflects the discrepancy already noticed for the EOS and 
mainly for the symmetry energy. This point will be dis- 
cussed in more details in a forthcoming paper (Baldo et 
al., in preparation). 

Our calculated BHF configurations with three-body 
forces are consistent with the measured masses of X-ray 
pulsars and radio pulsars (van Kerkwijk et al. 1995). 



4. Conclusions 

In conclusion, we computed some properties of NS's on 
the basis of a microscopic EOS obtained in the frame- 
work of BBG many body theory with two and three 
body nuclear interactions. Our Brueckner-Hartree-Fock 
EOS with three-body forces satisfies the general physical 
requirements (points i iv) discussed in the introduction. 
This is the main feature which distinguishes our EOS's 
with respect to other microscopic non-relativistic EOS 
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Fig. 7. The gravitational mass Mq, expressed in units of the solar mass A/q, is displayed vs. radius R (panc;l (a)) and 
the central number density ric (panel (b)). The notation is the same as in Fig. 3. The squares represent the values of 
the limiting configuration. 



(Wiringa et al. 1988, Engvik et al. 1994, 1996). Two dif- 
ferent two-body potentials have been used, the Argonnc 
AV14 and the Paris interactions. The latter turns out to 
produce a slightly stiffer EOS already at two-body level. 
The same three-body force, within the Urbana model, has 
been added to both two-body interactions, and this pro- 
duces similar saturation points close to the empirical one. 
The BHF calculations with three-body forces show again 
that the Paris potential gives an EOS stiffer than the one 
from the Argonne AVI 4 and a larger symmetry energy as 
a function of density. As a consequence the proton fraction 
is larger for the Paris and the Urea process onset appears 
at lower density. 

The calculated maximum mass is in agreement with 
observed NS masses (van Kerkwijk et al. 1995). However 
the limiting mass, radius and central density are relatively 
close to each other for the two interactions. In particu- 
lar the critical masses differ by no more than 8%, while 
the radius and central density are essentially equal. On 
the contrary all these three critical physical parameters 
are appreciably different from the ones derived from the 
Dirac-Brueckner approach (Li et al. 1992) as well as from 
the variational results. The latter displays also a rather 
different trend of the symmetry energy. 

For the Argonne AV14 (Paris) potential the BHF ap- 
proach with three-body forces predicts a lower limit for 
the mass above which neutron stars can support the di- 
rect Urea process, M^'^"" = IAMq (1.24M0). The stars 
with mass larger than this value cool very rapidly or not 
depending on the properties of nuclear superluidity (val- 
ues of the superfluid gaps, critical temperatures, density 



ranges for the superfluid transition) (Page et al. 1992, Baldo 

et al. 1995). The EOS's developed in this paper offer the 
possibility for a selfconsistent microscopic calculation for 
both the neutron star structure, and nuclear superfluid 
properties within the same many-body approach and with 
the same nuclear interaction. 
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